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Abstract 

We analyze the existence of localized waves in the vicinities of the interface between two dielectrics, provided one of them is 
uniaxial and lossy. We found two families of surface waves, one of them approaching the well-known Dyakonov surface waves 
(DSWs). In addition, a new family of wave fields exists which are tightly bound to the interface. Although its appearance is clearly 
associated with the dissipative character of the anisotropic material, the characteristic propagation length of such surface waves 
might surpasses the working wavelength by nearly two orders of magnitude. 
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1. Introduction 


Surface plasmons (SPs) appear in the disruption surface of 
isotropic media where the material permittivity changes of sign 
occurring with a metal in contact with a dielectric [1]. The rele¬ 
vance of these surface plasmons falls not only upon its inherent 
subwavelength localization, but also they enable an amplifica¬ 
tion of evanescent signals traveling near the surfaces |0]. These 
functionalities are being implemented within the lasts years for 
applications including optical sensing 0], signal filtering 0], 
subdiffraction focusing and subwavelength resolution imaging 

mum. 

Alternatively, it has been theoretically and experimentally 
demonstrated the existence of lossless surface waves at the in¬ 
terface of two different transparent dielectrics, provided one of 
them is anisotropic EMI. In opposition with SPs, this sort of 
surface waves have the peculiarity of possessing hybrid polar¬ 
ization [ 12]. The presence of hybrid surface waves with some 
parallel characteristics, additionally, may be found replacing 
the uniaxial medium by a biaxial crystal llj, 14], an indefinite 
medium Es, i6], and a structurally chiral material m, is]. 
The use of structured materials with extreme anisotropy offered 
an alternative to increase the range of directions of DSWs sub¬ 
stantially, as it is compared with the rather narrow range ob¬ 
served with natural birefringent materials [ 19, 20]. In partic¬ 
ular, striking results are attained if the anisotropic structures 
include metallic nanoelements, as it occurs for example with a 
simple metal-dielectric (MD) multilayer, a case where the angu¬ 
lar range may surpass half of a right angle I21.l22ll. Caused by 
the specific damping capacity of metals, however, the propaga¬ 
tion length of these DSWs is strongly limited by the penetration 
depth inside the lossy metamaterial El. 
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In this paper we perform a thorough analysis of DSWs tak¬ 
ing place in lossy uniaxial metamaterials. Special emphasis is 
put when the effective-medium approach (EMA) induces sat¬ 
isfactory results. We examine lossy metamaterials that exhibit 
closed spatial-dispersion curves, in the same manner that oc¬ 
curs with natural birefringent crystals. Contrarily, the intro¬ 
duction of losses leads to a transformation of the isofrequency 
curves, which deviates from spheres and ellipsoids, as com¬ 
monly considered by ordinary and extraordinary waves, respec¬ 
tively. As a consequence, two families of surface waves are 
found. One family of surface waves are directly related with the 
well-known solutions derived by Dyakonov [ 10]. Importantly, 
we reveal the existence of a new family of surface waves, which 
is closely connected to the presence of losses in the uniaxial ef¬ 
fective crystal. The dominant diffusion dynamics of the latter 
surfaces waves is thoroughly examined. 

2. Theory and configuration 

We first analyze the dispersion properties of DSWs supported 
by a ID array of semi-infinite layers made of a metal and a di¬ 
electric, displaced alternately as shown in Fig. [Ha), where w m 
(Wd ) is the width of the metallic (dielectric) layer within a given 
unit cell, and e m (6^) represents the dielectric constant of this 
material. In the semi-space z > 0 we have a lossless dielectric 
environment of relative permittivity e. An effective medium ap¬ 
proach based on the long-wavelength approximation is used to 
calculate the permittivity of the anisotropic metamaterial along 
its optical axis ey, that is the v-axis, and the permittivity in its 
normal direction, e ± . An schematic of the modeled uniaxial 
metamaterial is shown in Fig. Qlb). Under these conditions, 
DSWs of propagation constant kosw can be observed provided 
that e_L < £ < £|| Hot] . 

Form birefringence in the anisotropic metamaterial is es¬ 
timated in a simple way within the long-wavelength regime 
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Figure 1: (a) Schematic arrangement under study, consisting of a multilay¬ 

ered metal-dielectric metamaterial (z < 0) and an isotropic lossless medium 
(z > 0) of dielectric constant e. (b) The effective anisotropy of the stratified 
metamaterial is represented by the on-axis permittivity, e\\, and the permittivity 
over its perpendicular direction, e ± . A Dyakonov surface wave propagates in 
the vicinities of the interface between the isotropic and the effective uniaxial 
medium with propagation constant kpsw- 


that enables an homogenization of the structured metamate¬ 
rial [[24, [25]. The EMA demonstrates to be reliable for pho¬ 
tonic structures including elements with sizes that are signif¬ 
icantly smaller than the wavelength. At infrared and visible 
wavelengths, the skin depth of noble metals is clearly sub wave¬ 
length and, in this case, material homogenization requires that 
the metallic units had sizes of a few nanometers I2fj[27|,|28[]. In 
this case, the superlattice behaves as a uniaxial crystal whose 
optical axis is normal to the layers. The model estimates the 
relative permittivities along the optical axis, 


Cm 


" (1 -f)€ m +f€d' 

and transversally, 


( 1 ) 


c± = (l -f)c d + fe m , ( 2 ) 

where f = w m / iw d + w m ) is the metal filling factor. If we ne¬ 
glect losses by setting Im(6||) = 0 and Im(6 ± ) = 0, the effective 
birefringence of the MD superlattice is An = ^[e\\ - A 

small filling factor of the metallic composite may lead to an 
enormous birefringence, as inferred from Fig. [2] 

The effective permittivities of a MD periodic multilayer 
made of silicon dioxide ( e d = 2.25) and silver ( e m = -11.7 + 
/0.83 taken from [29]) at a wavelength To = 560 nm are dis¬ 
played in Fig. [2] showing how they vary for a different metal 
filling factor /. For the sake of clarity, we ignored the dissi¬ 
pative effects in the uniaxial metamaterial, and here we consid¬ 
ered the real part of e\\ and e± obtained from Eqs. ([]]) and ©. In 
this case, these are real and positive permittivities provided that 
/< 0.161, in addition leading to positive birefringence. Other¬ 
wise, the permittivity e ± became negative for a higher metal fill¬ 
ing factor. In addition, we considered cases where / > 0.0896, 
that is uniaxial metamaterials whose perpendicular permittivity 
satisfies e ± < e enabling the existence of DSWs; in this numer¬ 
ical simulation we considered e - 1 corresponding to air. 

In Fig. [3] we represent the solutions of the Dyakonov equa¬ 
tion, namely 


(k + Ke) (k + K 0 ) ( 6Ko + e ± K e ) = (fj| - e) (e - e x ) k%K„, (3) 



Figure 2: Effective permittivities e\\ and e± of a SiCF-Ag multilayered meta¬ 
material, at To = 560 nm, as a function of the metal filling factor /. In this 
simulation we disregard losses in the metamaterial, that is we used Re(e||) and 
Re(eD- The shadowed region indicates the range of values of / providing a 
necessary condition 0 < e± < e < 6y for the existence of DSWs, where we 
considered e = 1 corresponding to air. 



Figure 3: Isofrequency curve in the plane k x -k y of the solutions of the Dyakonov 
equation 0 calculated for e\\ = 2.625, 6 = 1, and e ± = 0.576, at a wave¬ 
length To = 650 nm. We include the dependence of k, k 0 , and K e on the spatial 
frequency k x . The spatial spectrum that corresponds to DSWs is bounded by 
k x = 0.619ko (K e = 0) and k x = 0.763ko (k = 0). 


providing the in-plane wave vector ( k x ,k y ) of the hybrid- 
polarized surface wave; note that the modulus of such a wave 
vector is the propagation constant kosw of the DSW. The elec¬ 
tromagnetic fields are evanescent in the isotropic medium, pro¬ 
portional to exp (-/c|z|), where ko = 2n/Ao. On the other side 
of the boundary, the ordinary and extraordinary waves in the 
effective uniaxial medium also decay exponentially with rates 
given by k q and K e , respectively. To illustrate the characteristic 
dispersion of DSWs, we solved Eq. ^ for the uniaxial Si02- 
Ag metamaterial described above, where the metal filling factor 
/ = 0.12, which is covered by air; again we assumed that the 
wavelength To = 650 nm. In this case, DSWs propagate within 
a wide angular region A ip = <£ max - <£ m j n , where tp is the po¬ 
lar angle of the surface wave vector ( k x , k y ). Specifically, we 
estimate A(p = 16.3° around a mean angle ip = 48.5° (DSWs 
lying in the first quadrant). The dispersion curve for DSWs is 
drawn in Fig. [3] As shown in that figure, <£ m i n is attained un¬ 
der the condition k = 0 (green solid line). On the other side 
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of the dispersion curve, <£> max is determined by K e = 0 (brown 
solid line) for which the extraordinary wave ceases to attenuate 

at v —>• +oo EMI. 


3. Simulation results and discussion 

Metal losses were disregarded in the analysis given above. 
In a realistic examination, however, the permittivity e m of the 
metal is not a real-valued constant, thus leading to a non¬ 
vanishing imaginary part of the permittivities e\\ and e ± of the 
metallic compound. The same occurs for the spatial frequen¬ 
cies k x and k y that are solutions of Eq. ©. The following pro¬ 
cedure can be used to find solutions of the Dyakonov equation: 
for a given real-valued spatial frequency k x along the optic-axis 
of the metamaterial, we calculated the complex-valued spatial 
frequencies k y that satisfies Eq. ©. This method is valid pro¬ 
vided that the excitation of the surface waves, given in terms of 
a plane-wave Fourier expansion, are produced at a plane paral¬ 
lel to y = 0, and that Im (k y ) > 0 [23]. As a result of utilizing 
a lossy uniaxial material, surface waves cannot propagate in¬ 
definitely and DSWs decay with a propagation length given by 

JUS] 

Ldsw = [2Im(^)] _1 . (4) 
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Figure 4: Dispersion curves for solutions of Dyakonov equation 0 considering 
a lossy metamaterial compound of permittivities e\\ = 2.625 + z'0.005 and e± = 
0.576 + z'0.100 at a wavelength Ao = 560 nm. We observe the existence of two 
families of hybrid-polarized surface waves, for which the real part of k, k () , and 
K e has a positive value, as shown in the insets. 

In Fig. 0] we represent the spatial frequency Re(fcy) of dissipa¬ 
tive DSWs propagating at different on-axis spatial frequencies 
k x over a lossy Si 02 -Ag metamaterial, where the metal filling 
factor / = 0.12. Focalization of the wave field in the vicinities 
of the plane z = 0 involves that the real part of k, k q , and K e has 
a positive value, as shown in the insets of Fig. 01 In principle, 
no restrictions must be imposed to their imaginary part Bill . 
Physically-realizable solutions of Eq. © are given within two 


spectral domains. The first domain has a spatial bandwidth that 
is comprised between k x = 0.615ko and k x = 0.753ko, where 
Re(ky) takes values from 0.639ko and 0.895ko, which is in good 
correspondence to the DSWs shown in Fig. 0] Interestingly, a 
second domain for Dyakonov solutions emerges within a higher 
bandwidth (of the spatial frequency k x ) ranging from 0.220ko to 
0.746ko. This new family of surface waves, here coined as type 
II DSWs, turns out to be the main result of our study. 
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Figure 5: Electric field \E y \, in arbitrary units, in the vicinities of the interface 
z = 0 for Dyakonov-type surface waves propagating at a spatial frequency: (a) 
k x = 0.67ko associated with the first family of DSWs, (b) k x = 0.67&o and (c) 
k x = 0.40ko corresponding the the new family of surface waves. 


Figure [5] depicts the modulus of the y-component of the elec¬ 
tric field, \E y \, in the vicinities of the interface z = 0 for hybrid- 
polarized surface waves propagating at a spatial frequency: (a) 
k x = 0.67 (k y = 0.752 + /0.0847 expressed in units of ko) asso¬ 
ciated with (type I) DSWs, (b) k x = 0.67 (k y = 0.397 + *0.136) 
and (c) k x = 0.40 (k y = 0.562 + /0.035) corresponding the the 
new family of surface waves (type II DSWs) that we found in 
our study. The x- and z-component of the electric field (not 
shown in the figure) are also localized around (z = 0). The 
analytical expression of the electric field for type II DSWs is 
the same as the expression that describes type I DSWs, which 
can be found for instance in Refs. [[l2] and JH. We observe 
a strong field confinement near the isotropic-anisotropic inter¬ 
face in all cases, though type I DSWs present a superior per¬ 
formance. This fact can also be inferred from the higher val¬ 
ues of the evanescent decay rates calculated for type I and type 
II DSWs, which are shown in the insets of Fig. 0] In addi¬ 
tion, a modulation on the tail of the evanescent field inside the 
lossy anisotropic medium is evident. This is caused by the ex¬ 
istence of a non-vanishing imaginary part of the decay rates k 0 
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and K e \ on the contrary such a modulation cannot be observed 
in the isotropic medium since field penetration is performed at 
a single rate k. Finally, the modulation of the field intensity is 
proportional to cos[Im(/c 0 - K e )x\, which can be observed pro¬ 
vided that its period is shorter than the penetration depth in the 
medium, as shown in Figs. Ob) and (c). 



Figure 6: Propagation length of type I DSWs and type II DSWs. 

In Fig. [6] we show the propagation length, defined in Eq. ©, 
and evaluated for type I DSWs and type II DSWs. Contrarily to 
what one may know by intuition, type II DSWs can propagate 
for a longer distance than type I DSWs, occurring in the spec¬ 
tral interval comprised between k x = 0.62ko and ko = 0.65ko. 
Moreover, the propagation length of type II DSWs can even 
surpasses 15 /im at low values of the spatial frequency k x . This 
is caused by the fact that the real part of k is very low within 
this spectral region, leading to low-decaying evanescent tails in 
the isotropic medium; see Fig. 0c) as an illustration of such an 
effect. Although the peak is placed close to the surface z = 0, 
in this case the average location of the field intensity is severely 
shifted toward the lossless material. Note that this also happens 
in long-range surface plasmons polaritons | ^]. 

4. Conclusions 

We numerically investigated the dispersion properties of 
Dyakonov surface waves propagating on the interface between 
an isotropic loss-free medium and a lossy uniaxial crystal, the 
latter resulting by applying the long-wavelength approach to 
metal-dielectric multilayered metamaterials. We conclude that 
the presence of absorption in the anisotropic medium limits 
the propagation length of these surface waves, but the spa¬ 
tial dispersion and field distribution near the material discon¬ 
tinuity varies slightly. In addition, we found a new family of 
electromagnetic waves which are confined near the isotropic- 
anisotropic interface and propagate obliquely to the optic axis 
of the birefringent material. These surface waves, here coined 
as type II Dyakonov surface waves, cannot be found in a config¬ 
uration using lossless materials. Unexpectedly, the propagation 
length of such surface wave fields might surpasses that related 
to type I DSWs. Potential applications of type II DSWs include 
ultra-compact waveguiding and biosensing. 
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